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Two out of the 18 structure functions describing 
the semi-inclusive lepton-nucleon deep-inelastic scatter- 
ing (SIDIS) process P, 0, see Fig. [U are associated at 
leading order of the hard scale Q with naively time- 
reversal odd (T-odd) transverse-momentum dependent 
parton distributions (TMDs), the Sivers and Boer- 
Mulders [1] functions. Their existence is ultimately re- 
lated to initial and final state interactions in QCD Q 
encoded in appropriately defined Wilson lines Hi • T- 
odd TMDs have unusual "universality" properties, and 
are predicted to have opposite signs Q in SIDIS and 
the Drcll-Yan process. The basis for this description is 
a generalized factorization approach which applies when 
the final state transverse momentum is small compared 
to the hard scale H^H, i.e. Ph± < Q in SIDIS. 

Data on the Sivers Boer-Mulders effect from SIDIS are 
available or forthcoming (Tl| - |2^ (for a recent review, see 
Ref. d^). Both effects were subject to p henomenological 
studies in SIDIS H-HI and DY Hiill . General aspects 
of the Sivers and Boer-Mulders functions were discussed 
i3-|4 
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.and quark model calculations were reported 
53I . Common to the quark model approaches is 



that the final (in SIDIS) or initial (in DY) state interac- 
tions are modeled by means of a one-gluon exchange (in 
(5^ steps were made to go beyond that). 

Among the most recent studies are the calculations in 
the light-cone constituent quark model (LCCQM) [53| . 
On the basis of the onc-gluon-cxchangc approximation 
formulated in the light-cone quantization formalism (s^ 
the Sivers and Boer-Mulders were modeled in [HI] us- 
ing light-cone wave-functions which were previously ap- 
plied with success to calculations of transversity, elec- 
troweak properties of the nucleon, generalized parton dis- 
tributions, distribution amplitudes [55l - l57| . and T-even 
leading- twist TMDs The results illustrate the rele- 
vance of different orbital angular momentum components 
of the nucleon light-cone wave function. 
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FIG. 1: Kinematics of the SIDIS process IN I'hX. The 
azimuthal angles of produced hadron and nucleon polariza- 
tion vector are <j>h and (ps- The transverse momentum of the 



hadron is Ph± ^ Q where Q 



{l-l'f 



The T-cven TMDs from this model (58| were shown 
to yield results in satisfactory agreement with available 
SIDIS data [5^. The purpose of this work is to explore 
whether also the T-odd TMDs from this approach [s^ 
can explain the corresponding SIDIS data. It is impor- 
tant to consider that in contrast to previous studies |55l - 
[ssj , we here probe more than the modeled nucleon light- 
cone wave- functions. In the present work we moreover 
also probe to which extent the "one-gluon-exchange ap- 
proximation" for T-odd TMDs works. 

The article is organized as follows. In Sec.|lT]the SIDIS 
process, and relevant obscrvablcs are briefly presented. 
In Sec. Uni the calculation of the T-odd TMDs [H is 
reviewed, positivity and evolution effects are discussed, 
and limitations of the approach conservatively disclosed. 
In Sccs. lIVI andlVlwe discuss the numerical results for the 
asymmetries and compare them to SIDIS data. Conclu- 
sions are drawn in Sec. IVII The appendix contains some 
remarks about the corrections due to the Cahn effect in 
the cos(20;i)-asymmetry in unpolarized SIDIS. 
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II. SIVERS AND BOER-MULDERS IN SIDIS 

The kinematics of the SIDIS process, the momenta 
I, l\ q, Q2 and Phi_ = \Ph±\ are defined in Fig. [T] The 
SIDIS variables aie x = /{2P ■ q), y = {P ■ q)/{P ■ I), 
and z = {P ■ Pfi) / {P ■ q) where P is the nucleon momen- 
tum. Denoting by (Tq its spin- and ^/j-independent part, 
the SIDIS cross-section a (differential in x, y, z and the 
azimuthal angle (f>h which we do not indicate for brevity) 
can be written as 



dV = d^ 



sin 
UT 



(1) 



where Pi (y) = (1 — — y + ^ 2/^) up to (systematically 
neglected) power-suppressed terms, and the dots indicate 
terms due to T-even TMDs or subleading- twist [60| . 



jcos(2(/)h) 
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/F(ju and similarly Ajjj: 



i<Ph-<t>s) 



are the azimuthal asymmetries defined in terms of struc- 
ture functions. Fm/ is the unpolarized structure func- 
tion which gives rise to ctq- Here the first index U de- 
notes the unpolarized electron beam. The second index 
U/T refers to the target polarization which can be un- 
polarizcd/transvcrsc with respect to the virtual photon 
(in experiment the transverse target polarization is with 
respect to the beam, of course, but this is up to power- 
corrections the same). The superscript reminds us of the 
kind of angular distribution of the produced hadrons with 
no index describing an isotropic 0?i-distribution. In the 
Bjorken-limit the relevant structure functions are given 
at tree-level by [i^l 
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(2) 
(3) 
,(4) 



where h = Pfi±/Ph± and M {rah) is the mass of nucleon 
(produced hadron). The convolutions arc defined as 



C 



w i J 



X w 



{Pt, KT)Y.elxr{x,pT) r{z,KT), (5) 



with a generic TMD and transverse momentum de- 
pendent fragmentation function J", and pt = |PtI and 
Kt = \Kt\- In Eqs. pH U) D1 is the unpolarized and 
H^'^ the Collins [6l|, [63 fragmentation function. Notice 
that in this tree-level treatment one neglects soft factors 



In order to solve the convolution integrals we will make 
use of the Gaussian Ansatz. This step could, in principle, 
be avoided by considering adequately weighted asymme- 
tries [1] . It is given by 



3''{x,pT)=f{x) 



exp(-pj/(p|(j))) 



r(z,Kj 



Independently of the model for transverse momenta 
Fuu{x, z) = Yl,a ^ /i (^) -^1 (^)i whilc with the Ansatz 
(|6l) we obtain for the structure functions [531 



F-^^^-''^\x^z) = -B,Y.^lxfti'^%x)Dt{z)^ (7) 
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{x,z) = B,Y.eixh^^'^\x)Ht^"'^'\z), (8) 
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Bo = 



z 0rM 



{z^l^UftT)) + {KUD,))Y/^-' 
8zM[tt{K^{H^))]-^/^ 



(9) 
(10) 



l + z^{p^hi))/{KUH,))' 
with the transverse moments defined as 

2A/2 

J^(i/2)-(z)= f d^K^^J^^z^KT). (11) 

J ZZTflji 

The Gaussian Ansatz (jH) is for /f and Dt phenomeno- 
logically well supported for {Pha.) ^ Q [SJ, 113]. For 
polarized TMDs it is supported approximately by some 
models [H, [H, [6^ 111]. A rigorous description of pt- 
effects would require methods along the QCD-bascd for- 
malism of 67 1 , as implemented in (gsI . |69| , but for our 
purposes the Ansatz ([6]) will be suflicient. 

It is important to stress that the description of the 
SIDIS process in the Bjorken-limit within TMD factor- 
ization is generically valid only up to p ower corrections 
suppressed by the large scale Q p, [lOl • Typically such 
corrections do not factorize (TMD factorization beyond 
leading twist is presently unclear [t^I in SIDIS) and must 
be excluded experimentally by studying the Q-behavior 

icos 20/1 



of the observables. In the case of Fjjjj ^'^ the parton 
model provides a way to estimate one of the possible 
power corrections which is due to the Cahn-effect [7l| . 
This effect gives rise to cos0/j [13: 113 and cos{2cj)h) mod- 
ulations in the unpolarized SIDIS cross section which are 
suppressed by respectively {pT{fi))/Q and (Pt(/i))/Q^- 
In particular, the latter contributes a power-correction to 
the Boer-Mulders asymmetry, which is not n egligibl e in 
the kinematics of present SIDIS experiments [18|, |l3 • 

The Cahn-effect and its power correction to A'^^^'^'^\ 
however, are to a good approximation fiavor independent 
[U, 113]. Below in Sec. |V]we will therefore consider the 
difference of the asymmetries for tt~ and 7r+. Further 
details are discussed in the Appendix. 
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III. T-ODD TMDS IN THE LCCQM 

In this section we briefly review the calculation of the 
T-odd TMDs in the light-cone constituent quark model 
[53j . Within this model, the gauge-link operator entering 
the quark correlation function which defines the TMDs is 
approximated by taking into account one-gluon exchange 
between the struck quark and the spectator quarks in 
the nucleon described by (real) light-cone wave functions 
(LCWFs). Working in the light-cone gauge = 0, the 
T-odd quark TMDs are written in terms of overlap of 
LCWFs, convoluted with the gluon propagator obtained 
from the expansion of the gauge link. The final result 
is of order as, which enters as an overall multiplicative 
factor and is understood at the initial scale fiQ of the 
model. 

The value of as{no) was fixed by determining the ini- 
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tial model scale /Xg following [72| as follows. Since the 
model has only valence quark degrees of freedom, one 
requires to be such that evolving the (total in the 
model) momentum fraction carried by valence quarks, 
{x{fJ.o))v = J2ql da;a;(/f - ff){x,fio) = 1, from fiQ to 
experimentally relevant scales one matches the NLO phe- 
nomenological value {x{Q))v = 0.36 at = lOGeV^ 
[73} . In previous works fj-o was tuned to reproduce exactly 
the phenomenological value of {x{Q))v- Here we content 
ourselves with reproducing it within 10%, which is ac- 
ceptable in view of the generic model accuracy of (10- 
30) % ii]. Allowing for a 10 %-overestimate of {x{Q))v 
yields a higher ^0 and smaller as{no), and better con- 
vergence. The evolution was performed in the variable 
fiavor-number scheme with heavy-quark mass thresholds 
rric = 1.4 GeV, m& = 4.75 GeV, = 175 GeV with 
a^L°(M|) = 0.12018 at next-to-leading order (NLO, MS 
scheme) from [t^. We obtain for the initial model scale 



^NLO ^ 508 MeV, 



a 



NLO 



47r 



0.128 



(12) 



(notice that A^^lo^^ = 402, 341, 239 MeV in [zl). The 
strategy of fixing the model scale is basically the same as 
in previous works. Here we introduced the variable flavor 
number scheme, and updated the value of as{M'^). As 
a result, we find a higher value for the hadronic scale of 
the model, corresponding to a somewhat smaller coupling 
compared to [H^l- For the calculations of T-odd TMDs, 
we will adopt the value of as in Eq. ([T^ . 

The LCWF of the nucleon is modeled as described in 
[53I [sst. In particular, to disentangle the spin-spin and 
spin-orbit quark correlations encoded in the TMDs, we 
expand the three-quark LCWF in a basis of eigenstates 
of orbital angular momentum, which yields 6 indepen- 
dent amplitudes corresponding to different combinations 
of quark helicity and orbital angular momentum [t^ . 

Assuming SU(6) symmetry, these light-cone ampli- 
tudes have a particularly simple structure, with spin 
and isospin dependence factorized from a momentum- 
dependent function with spherical symmetry. Under this 
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FIG. 2: Transverse moments Svr^'^^'^) ^'^'^ h^'^^^'^ix) in a 
proton for up (left) and down (right) quark, as function of x. 
The dashed curves show the results at the hadronic scale /^o 
of the model . The solid curves correspond to the results after 
LO evolution to = 2.5 GeV'^, using the evolution patterns 
of the unpolarized parton distribution (transversity) for the 
Sivers (Boer-Mulders) function. 



assumption the orbital angular momentum content of the 
wave function is fully generated by the Melosh rotations 
which boost the rest-frame (canonical) spin of the quarks 
into light-cone helicities. For the momentum-dependent 
part of the LCWF we adopt the phenomenological de- 
scription with parameters fitted to hadronic structure 
constants fro m 17611 which gave satisfactory results in pre- 
vious works |56l459| . 

The results for the (l)-momcnts of the Sivers and Boer- 
Mulders functions as defined in pT|) are shown in Fig. [5] 
The dashed curves show the results at the hadronic scale 
of the model, while the solid curves are obtained by ap- 
plying LO evolution to = 2.5 GeV^ Since the 
exact evolution equations for the T-odd TMDs arc still 
under study [TSl - fsll . [sst . we use those evolution equations 
which seem most promising to "simulate" the correct evo- 
lution. 

The chiral-even Sivers function is evolved according to 
the LO non-singlet evolution pattern of /f (x), which has 
the advantage of preserving the Burkardt sum rule . 
This non-trivial constraint for model calculations is sat- 
isfied in the LCCQM at the initial model scale [s^l- For 
the chiral-odd Boer-Mulders function we use the evolu- 
tion pattern of h1{x). The so evolved model results are 
consistent with the first extractions [26l - l3l| concerning 
signs and magnitudes of the various fiavors, and the po- 
sitions of their maxima in x. 

As can be seen in Fig. [2l evolution effects are impor- 
tant, and this raises the question how reliably we esti- 
mated them here. The exact answer to this question will 
be given after the TMD-evolution formalism discussed 
in [63] for fi{x,pT) has been generahzed to include the 
Sivers and Boer-Mulders effects. Until then, we can try 
to gain some intuition by investigating how much the re- 
sults would change if we used different ways of simulat- 
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ing the evolution effects for T-odd TMDs. For tfiat we 
explored also alternative evolution patterns, and found 
that the corresponding results vary within (10-20) %. 
This indicates that the theoretical uncertainty associated 
with evolution effects is not larger than the model accu- 
racy [i^l ■ An interesting result is that the model supports 
approximately the Gaussian Ansatz, both for T-odd [1^ 
and T-even [s!] TMDs. 

A point deserving particular attention are inequalities 
among twist-2 TMDs [i^, which can be written as 

^Siv(2^'-PT) = fl{x,pTf - gl{x,pTf 

9t^{^,PTf - ^ ftr^^^PT? > , (13) 

-^hti{x,PTr-^ht\x,PTr>o. (14) 



In a consistent model framework the inequalities (|13p- 
([13]) should hold. However, as observed in [s^, for cer- 
tain values of x and px numerous models violate already 
the weaker inequalities \ f^{x,pT)\ < fi{x,PT), and 

^ \ht''{x,pT)\ < f!ix,PT) following from (USD, (HI. 

The reason is apparent. The calculations of T-odd 
TMDs arc within the given models correct. But while 
T-odd TMDs arc evaluated to 0{as), the "expansion" of 
T-even TMDs is truncated in the quark models at 0{a^). 
In order to preserve unitarity, and hence the inequalities, 
one should evaluate also T-even TMDs consistently to 
order 0{as). To best of our knowledge, this has so far 
not been done in any model, and would also go beyond 
the scope of the present work. 

Having established that one cannot expect the inequal- 
ities (|T3lfT4)) to be satisfied in our approach either, let us 
see for which values of x and pr they are violated. For 
that, in Fig. [3] we plot 7'|i^,BM(2^'PT), Eqs. dHIHl), vs. 
p^ for selected values of x = 0.2, 0.3, 0.4 at the low model 
scale. Vg^^ bmC-^'Pt) should be always positive, and we 
see that this condition is violated only at small x and for 
p^ significantly larger than the respective {px)- Now it is 
important to recall several facts. First, quark models do 
not describe reliably small a;. Notice that the region of 
x < 0.2 at the low model scale corresponds after (correct 
or approximate) evolution to x < O(10~^). Thus, prac- 
tical consequences (if any) of the violation of inequalities 
are beyond the region of x where quark models can be ap- 
plied. Second, in the following we will restrict ourselves 
to the use of pr-integrated results for the (l)-moments, 
and those are dominated by the regions of pT of the order 
of magnitude of (p^) and smaller, i.e. where the inequal- 
ities are satisfied. 

To conclude, though we have to admit that the model 
predictions do not satisfy the inequalities for all x and pr , 
we can be assured that — in the way we will use them 
for phenomcnological calculations — this will not yield 
any artifacts. In this sense we will consider our results 
as compliant with positivity constraints. 
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FIG. 3: The positivity relations involving the Sivers (up- 
per panels) and Boer-Mulders functions (lower panels) from 
Eqs. (|13p and (|14|l . respectively, as function of at different 
values of a;: x = 0.3 (short-dashed curves), x — 0.4 (long- 
dashed curves), and x = 0.5 (solid curves). The left (right) 
panels show the results for up (down) quark. 



IV. SIVERS ASYMMETRY 



In this Section wc discuss the asymmetry A^!^^'^''"''^^'' 
due to the Sivers effect. With the model results referring 
to a low scale, it is not only necessary to evolve the Sivers 
function in x, but also to account for pT-broadening (67j 
(all this applies equally to the Boer-Mulders function). 
In principle, one could feed the CoUins-Soper-Sterman 
(CSS) equation [gJ with the model results as initial con- 
ditions. We refrain from this step, because the CSS for- 
malism is not yet developed for cases including polariza- 
tion effects, while in the unpolarized case the result is 
known: one would expect to reproduce, within model ac- 
curacy, the phenomcnological value (Pxifi)) ~ 0.4 GeV^ 
at HERMES energies. Therefore we proceed as follows. 

We assume that the Gaussian shape is approximately 
preserved through the CSS evolution but the Gaussian 
widths increase with energy [g^, |6^ , and we assume that 
to lowest order approximation this pT-broadening is po- 
larization independent. Wc use the model predictions 
for f^r^^^'^(x) (approximately evolved in x as discussed 



in Sec. IIII[) which arc presumably less affected by Su- 
dakov effects [s^l than f^r^{x) with which one also could 
work within the Gauss model. In the expression ([7]) for 
the Sivers asymmetry we use {K^{Di)) = 0.16 GeV^ [l^l 
and the Gaussian width of the Sivers function, which 
is (pK/iV)) « 0.9(p|(/i)) in the model [H. Accord- 
ing to our assumption one may expect this prediction to 
be roughly valid also at experimentally relevant scales, 
because to lowest order approximation pT-broadcning ef- 
fects are polarization-independent. In [59[ positive expe- 
rience was made with such estimates of jir-broadening 
effects. Further studies are required for more precision. 
In the numerator of the Sivers asymmetry we use 
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FIG. 4: The single-spin asymmetry A^y"^'' '^'^^ for pion production oflF proton in SIDIS, as function of x. The HERMES data 
are from [l^. The theoretical curves are obtained on the basis of the LCCQM predictions for fyr^''''{x) [53^ . 



Af^VI-s' (a) Af^V*s' (b) Af^VI-s' (c) Af^^'I's' (d) 
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FIG. 5: yl'j^y''^'^"'^^^ for charged pion (hadron) production off deuteron (proton) in SIDIS, as function of x. The COMPASS 
data are from [T^, [ij . The theoretical curves are obtained on the basis of the LCCQM predictions for f^r^^'^ix) The 
dotted curve in panel (c) shows the results without the effects of pT-broadening (see text). 
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FIG. 6: ylj}'^^^" for pion production off neutron in SIDIS, as function of x. The theoretical predictions are obtained on the 
basis of results for fiT^'''^ix) from LCCQM for the kinematics of the Hall-A experiment at Jefferson Lab [2^ . 



model predictions for ff{x) LO-evolvcd to a scale of 
2.5 GeV^. For D1{z) we use the LO parameterization 
from [s^ at the same scale. 

In this way we obtain in the kinematics of the HER- 
MES experiment, (Q^) ~ 2.5 GeV^ and 0.2 < z < 0.7, 
the results for the x-dcpendence of the Sivers asymmetry 
for pion production off a proton target shown in Fig. |4l 
Keeping in mind that the quark model approach is not 
expected to be reliable in the small- a; region, we observe 
a good description of the data within the accuracy of our 
approach. (Due to the absence of sea quarks in our ap- 
proach the Sivers asymmetries for kaons would be very 
similar to the pion asymmetries. The explanation of a 
possible difference in tt"*" and Sivers asymmetries 15 1 



is beyond the scope of our approach.) 

Next we discuss the COMPASS data which have (Q^) 
similar to HERMES and 0.2 < z < 1, but were mea- 
sured at a significantly higher s = (P + /)^ w 2MAri?bcam 
with i^boam = 160 GeV (vs. Sbcam = 27 GcV at HER- 
MES). From Drell-Yan it is known that at higher en- 
ergies the TMDs tend to be broader Q as follows 
from the CSS formalism [13, . There are indications 
that this is also the case in SIDIS l641, namely one has 
{Phi) = 0.27 GeV^ at HERMES [8| to be compared 
with (P2^) ^ 0.41 GeV^ at COMPASS [H. Both re- 
sults refer to a common z ~ 0.55 but different {x). If we 
assume an x-indcpendcnt Gaussian width of /f (the data 
are compatible with this assumption the observed 
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FIG. 7: The difference of azimuthal asymmetries A'^^ for negative and positive pions or hadrons, as function of x. The 
experimental points were obtained taking the differences of preliminary tv~ and tt"^ HERMES data [1^, and preliminary h~ 
and h'^ COMPASS data [l^. The error bars show the propagation of statistical errors, and do not include systematic errors. 
The theoretical curves are obtained using hf^^^''{x) from the LCCQM [s^l- Panel (d) shows a prediction for Jefferson Lab. 



broadening of (P^^) ~ {K^{Di)) + has to be 

attributed to broadenings of the Gaussian widths of /f 
and Df. In order to crudely estimate these effects we 
assume that at COMPASS the widths of /f and are 
equally broadened compared to HERMES as follows 

(Pt(/i)) COMPASS = (Pt(/i)) HERMES + S{kx) , 
(A'f(£)l)) COMPASS = (-ft'|(/l))HERMES + S{k'^) , (15) 

with w 0.11 GeV^ needed to explain the larger 

{P^±) at COMPASS. This change in parameters also af- 
fects (broadens) the Gaussian width of the Sivers func- 
tion (estimated here as {Prifrr)) ~ 0-9 (PtI./!)): see 
above). In this way we obtain the results shown in Fig.[S] 
We observe a good agreement of the model results 



with the COMPASS data on the Sivers effect in tt^ 



pro- 



duction from a deuteron target [l3| in Figs. [5^, b, and 
charged hadron production from a proton target |14l | in 
Figs. [SJ:, d. (For simplicity we approximated the results 
for charged hadrons by charged pions, which account for 
about 90% of charged hadrons at COMPASS energies.) 

Notice that, had we neglected the p^-broadening at 
COMPASS as compared to HERMES, the model would 
have clearly overestimated 7r+ data from proton at COM- 
PASS, see Fig. [St. Indeed, Fig. [S}; seems to indicate that 
the p^-broadening could even be somewhat stronger than 
estimated on the basis of (|T5]) . In Figs. [5^, b, d the ne- 
glect of broadening effects would be less significant, or 
completely within error bars. 

Finally, we present predictions for a forthcoming Hall- 
A experiment at Jefferson Lab [2lj with a 5.9 GeV beam. 
Here s is sufficiently close to HERMES, such that it is 
not necessary to account for p^-broadening effects, in 
contrast to COMPASS, Eq. ^ [gll. In this experiment 
(g^) - 2 GeV^ and 0.42 < z< 0.66. The rcsuhs for this 
kinematics are shown in Fig. [51 At Hall- A the asymmetry 
is larger because roughly Bq{z) oc z in ([7]) and Hall- A 
probes larger z compared to HERMES and COMPASS. 

To summarize, we find that our model framework pro- 
vides a satisfactory description of the SIDIS data on the 
Sivers effect from HERMES and COMPASS [iMl. Fu- 
ture data from Jefferson Lab will allow further tests. 



V. BOER-MULDERS ASYMMETRY 

In this section we focus on the asymmetry A'^^^'^'^^^ 
due to the Boer-Mulders effect. The asymmetry ([5]) is 
calculated similarly to the Sivers asymmetry in Sec. IIVI 

We use the model predictions for h^'^'^^'^{x) from (ssj . 
For the Collins function and its Gaussian width we use 
the information from [s^, [s^] , and for the width of the 
Boer-Mulders function we use the LCCQM model pre- 
diction {Pxih't)) ~ 0-95 (pf.(/i)) nil, which we again 
assume to be approximately valid at experimentally 
relevant scales. For /f(x) in the numerator of the 
asymmetr y w e use the model predictions LO-evolved to 
2.5 GeV^ [9^, and for D^{z) the LO parameterization 
from [s^ at the same scale. The pr-broadening effects 
at the higher energies in the COMPASS as compared 
to HERMES, are estimated similarly to Sec. IIVI i.e. we 
use the HERMES value {K^{H^)) « Q.2b{K^{Dx)) from 



j87l |. and consider py-broadening analog to Eq. p^ . 

As discussed in detail in Sec. [Hi the Cahn-effect f7l| 
generates an l/Q"^ power-correction to the cos(20/i) mod- 
ulation in the unpolarized SIDIS cross section which can- 
not be neglected in the kinematics of the HERMES or 
COMPASS experiments [Ull^]. In principle, one can try 
to model this power-correction [3l|. For that one could 
use the updated phenomenological results for {Pxifi)) 
and {K^{Di)) [64 1 which are sufHcient to determine the 
Cahn-effect. Alternatively, one could explore data on 
the cos(2(/)/j)-asymmetry of neutral pions, where due to 
the Savor-dependence of the Collins function [H, |83| the 
leading-twist Boer-Mulders effect largely cancels, see Ap- 
pendix. In any case, this step constitutes an additional 
modeling step (of a non-factorizable twist-4 contribu- 
tion), and we prefer to avoid it. 



However, the Cahn effect contamination in 



IS 



flavor-independent to a good approximation, and largely 
cancels out in differences of 7r~ and 7r+ asymmetries. 
Such differences can be determined from the preliminary 
data [HjIi^ and we shall confront our model results with 
them instead. The results are shown in Figs. [7^-c. (We 
again approximate at COMPASS results by tt^-) 
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Again it is important to recall that the quark results 
should not be expected to be reliable in the small-a; re- 
gion. In the region of a; ^ 0.1 the model describes well 
the size of the asymmetry differences, but it seems not 
to follow the trend of the data at larger x > 0.2. How- 
ever, one has to keep in mind the preliminary status of 
the data [l^, [l^. Moreover, considering systematic er- 
rors of the data which are not included in Figs. [7^-c, the 
discrepancy could be well within model accuracy. 

From the available Jefferson Lab data [l^, [l^ the dif- 
ference of TT^ and 7r+ cos(2(/)/j)-asymmetries cannot be 
accessed, but forthcoming CLAS data will allow that [20| . 
In the CLAS experiment we have (Q^) ~ 1.9 GeV^ with 
0.1 < x < 0.6 and 0.4 < z < 0.7. The predictions for 
this kinematics are shown in Fig. [7Jl. Since the x values 
of the CLAS experiment cover the region dominated by 
the valence-quark contribution, these data will provide 
an important test of the model. 

To summarize, our approach is compatible with pre- 
hminary data from COMPASS and HERMES. Further 
insights can be expected from Jefferson Lab, before [l^l 
and after the 12 GeV beam energy upgrade [2^, and on 
long term from the future Electron- Ion Collider [2^ . 



VI. CONCLUSIONS 

In this work we have studied the leading-twist az- 
imuthal asymmetries in SIDIS due to T-odd TMDs on 
the basis of predictions from the light-cone constituent- 
quark model [53| . Since the model results refer at a low 
hadronic scale, we discussed how to take into account the 
effects of the evolution for the description of data refer- 
ring to high scales of typically several GeV^ . We tackled 
this issue in two steps. First, for the px-depcndence of 
the distributions we employed the Gaussian Ansatz and 
expressed the asymmetries in terms of (l)-transverse mo- 
ments of TMDs. The Gaussian widths of the distribu- 
tions were assumed x- independent, which is supported 
by phenomenology, and p^-broadening effects were es- 
timated. In the second step, we evolved the transverse 
moments of TMDs to experimental scales by employing 
those evolution equations which seem most promising to 
simulate the correct evolution. For the Sivers distribu- 
tion we used the non-singlet evolution pattern of fi{x). 
This allows to preserve the Burkardt sum rule — valid at 
the initial scale of the model [s^ — also at higher scales. 
For the chiral-odd Boer-Mulders function we used the 
evolution pattern of transversity. 

We obtained a good description of the Sivers asymme- 
try, and satisfactory results for the Boer-Mulders asym- 
metry in comparison with available experimental data 
from' HERMES and COMPASS. In the case of the Boer- 
Mulders asymmetry we considered differences of asym- 
metries for Tr~ and tt"*" to avoid the modeling of twist-4 
power corrections ('Cahn effect'). Furthermore, we pre- 
sented model predictions for forthcoming experiments at 
Jefferson Lab, which will extend the available data far 



into the valence- a; region where the model is expected to 
work best, and provide an important test of its dynamics. 

Our results indicate that the use of the one-gluon- 
exchange-mechanism to model T-odd TMDs (as imple- 
mented in [5^ ) yields phenemenologically reasonable re- 
sults, although a truncation of the expansion of the 
gauge-link at 0{as) seems not a priori justifiable. 

The present work completes the study in Ref. [H^ 
where leading-twist spin asymmetries due to T-even 
TMDs were calculated. We observe that the light-cone 
constituent-quark model, based on overlap representa- 
tion of TMDs in terms of light-cone wave functions, pro- 
vides a good description of intrinsic transverse parton 
momentum effects in the range of applicability of the 
model. 
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Appendix A: Remark on Cahn effect 

Finally, a remark is in order on A™^^^'^''^ for definite pi- 
ous, where the Cahn effect l/Q^-correction cannot be ne- 
glected Ull, [13] though it could be small in largest x-bins 
at COMPASS. In principle, one may model this contri- 
bution [m with due care to the energy dependence of the 
Gauss parameters [13, • However the problem remains 
how to test independently such an additional modeling 
of a presumably non-factorizing twist-4 term. We are not 
aware of a 'rigorous procedure', but the following obser- 
vation may turn out helpful. The favored (u — tt"*") and 
unfavored {u — >■ tt") Collins fragmentation functions are 
similar in magnitude but have opposite signs (sgI [STj . 

In the A^^^'^'*^"*^'' asymmetry for 7r° (which is due to 
the Collins effect and potentially not or less affected by 
power-corrections) this yields nearly exact fiavor cancel- 
lations, as seen in data (9l| . Since this is a property of 
the Collins function, one expects the Collins effect to also 
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cancel out in the n'^ cos(2(/)/i)-asymmetry, i.e. 

.COs(20h)/ 0\ ^ ACOs{2cl>h) 

^I7t/ I ~ ^UU, Calm ^ " ' 



(Al) 



where ■^uu'^C^n ^-'^^ Cahn effect contribution, which 



is negative [71[ and largely flavor independent [31|, l6J] . 
This would then mean that 

^C/j/ 17!^ J ~ ^UU \^ JBM + Cahn ' l^^J 

^cos(20h)/ ^ ACOs(24,h) r -\ , -cos(20fe) . . o^ 



This is not to be confused with relations due to isospin 



invariance, which allow one to express tt" cross sections in 
terms of tt^ cross sections (i.e. which connect the numera- 
tors or the denominators of the asymmetries). Eqs. (jAip - 
(|A3[) indicate how to model the Cahn-effect in a given 
experiment. This might be a more reliable procedure 
than using other sources of information. The procedure 
can be iteratively improved to take into account flavor- 
dependencies in the Cahn-efFect and non-exact cancella- 
tions of the Collins effect in the neutral pion cos(2(/)/j)- 
asymmctry. In particular, Eqs. (|A2|) and (|A3|) show our 
underlying assumption in Sec. |V] that in the difference 
of charged pion cos(20/i)-asymmetries the Cahn effect 
largely cancels out. 
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